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An  expression  has  been  developed  for  the  longi¬ 
tudinal  component  of  the  vibratory  force  exerted  on  a  prolate 
spheroid  by  the  operation  of  a  marine  propeller  In  a  space- 
varying  field  (wake).  Two  evaluation  schemes  have  been  con¬ 
sidered:  one  by  Integration  of  the  pressure  signal  over  the 
surface  of  the  ellipsoid  and  the  other  by  means  of  Lagally's 
theorem  with  the  ellipsoid  represented  by  a  known  source- sink 
distribution.  Numerical  calculations  Indicate  the  Important 
role  played  by  propeller  clearance  and  slenderness  ratio  In 
the  magnitude  of  the  vibratory  force. 
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INTRODUCTION 


During  the  past  few  years,  theoretical  Investiga¬ 
tions  have  been  undertaken  at  Davidson  Laboratory  of 
Stevens  Institute  of  Technology  Into  the  vibratory  pressure 
and  velocity  field  around  an  operating  marine  propeller  with 
the  object  of  determining  the  forces  exerted  on  nearby 
bodies.  Early  studies^'®  were  restricted  to  the  case  of  uni¬ 
form  Inflow  to  the  propeller  which  Is  represented  by  a  llne- 
vortex  array  to  give  the  effect  of  loading  and  by  a  source- 
sink  distribution  for  the  blade  thickness  effect.  Later 
studies  of  the  propeller  field  treated  nonuniform  Inflow  con¬ 
ditions  since  In  the  case  of  a  marine  propeller  located  be¬ 
hind  a  hull,  nonuniformity  of  the  Inflow  field  Is  a  common 
occurrence.  In  these  studies®*'*  which  were  concerned  with 
both  acoustic  and  hydrodynamic  media,  the  propeller  was  rep¬ 
resented  by  axial  and  tangential  doublet  distributions  on  the 
propeller  blade  axis  with  strength  depending  on  radial  and 
angular  position.  In  considering  propeller  operation  under 
nonuniform  inflow  conditions  (wake),  the  analysis  takes  cog¬ 
nizance  of  a  situation  more  realistic  than  the  open  water 
condition,  since  the  effect  of  the  boundary  Is  taken  implicit¬ 
ly  into  account  through  the  wake  formation.  In  other  investi¬ 
gations,  Initiated  by  Davidson  Laboratory  toward  fulfillment 
of  the  long-range  objective  of  the  series,  expressions  have 
been  developed  for  the  vibratory  forces  eind  moments  produced 
by  a  marine  propeller  on  a  doubly- Infinite  rigid  plate®  and 
an  Infinitely  long  rigid  strip,®  Results  of  the  latter  study 
are  surprisingly  simple  auid  Indicate  that  a  vibratory  force 
of  considerable  magnitude  could  be  obtained  from  the  treinslent 
field  generated  by  an  operating  propeller. 

The  present  study  endeavors  to  determine  the  vibra¬ 
tory  forces  exerted  on  an  ellipsoid  of  revolution  by  a  marine 
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propeller  operating  in  a  specified  wake  behind  a  prolate 
spheroid. 

The  study  has  been  carried  out  under  Bureau  of  Ships 
Fundamental  Hydrodynamics  Research  Program  S-R009  01  01, 
Contract  Nonr  263(  l6)  eind  administered  by  David  Taylor  Model 
Basin* 


QUASI- STEADY  AND  UNSTEADY  VIBRATORY 
FORCE  ON  A  PROLATE  SPHEROID 

Vibratory  Force 

The  propeller  disc  Is  located  at  a  distance  f  from 
the  center  of  the  spheroid  whose  major  axis  coincides  with 
the  longitudinal  axis  of  the  propeller  ( see  sketch) .  The 
prolate  spheroid  Is  defined  in  terms  of  Its  major  axis  2a  and 

Its  eccentricity  e.  Two  car¬ 
tesian  coordinate  systems  are 
used  with  the  same  direction 
of  positive  axes  but  with  dif¬ 
ferent  origin.  The  axes 
located  at  the  propeller  center 
are  designated  by  x,  y,  z  and 
the  axes  with  origin  at  the 
center  of  the  ellipsoid  of 

revolution  by  x^,  y.  r,. 

o  o  o 

The  relation  between  the  cartesian  coordinates 
y^,  and  ellipsoidal  coordinates  ju,  f  are  given  by 


kju? 

K  C<  < 

k(  l-iU^) 

(  ?®-l)  cos  f 

-1^M<1 

k(l-|u)'^ 

( ^^-l)  sin  f 

0^f^2iT 

Sketch  No.  1 
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For  a  given  ellipsoid  of  known  eccentricity  e  and  major 
semi- axis «  a, 

“  p  k  =  ae,  M  =  ^ 

1/2 

The  radial  distance  r^  =  (x®  +  y®)  '  =  a(l-e^)  (l  - 


(1) 


and  the  x^-coslne  direction  of  the  unit  vector  tangent  to 
the  surface  C  =  constant  is  given  by 


The  propeller  located  behind  a  ship  operates  in  space-varying 
Inflow  conditions  so  that  its  rotating  blades  experience  a 
time-dependent  gust.  On  the  basis  of  this  physical  reason¬ 
ing  an  approximate  theory  has  been  developed"^  utilizing 
two-dimensional  unsteady  theory  in  a  stripwise  fashion.  It 
is  true ,  on  the  other  hand,  that  the  pressure  of  the  pro¬ 
peller  disturbs  the  fluid  field  around  the  body  so  that  the 
stream  lines  adjacent  to  the  body  are  distorted.  This  amounts 
to  distortion  of  the  original  shape  of  the  body.  The  effect 
on  the  boundary  of  the  presence  of  the  propeller,  the  so- 
called  "image  effect"  necessary  to  restore  the  boundary,  was 
rigorously  evaluated  in  the  case  of  the  doubly- infinite  rigid 
plate®  and  later  for  an  infinitely  long  rigid  cylinder.®'® 

In  the  latter  studies  it  was  found  Independently  that  the 
vibratory  force  exerted  on  the  cylinder  due  to  the  "imege 
effect"  is  identically  the  same  as  that  developed  directly 
by  the  propeller  action.  To  the  best  of  the  authors’  know¬ 
ledge,  there  is  not  at  present  a  method  of  constructing  the 
"image  potential"  for  the  case  of  an  ellipsoid  of  revolution 
in  the  presence  of  the  propeller.  This  case  is  complicated 
further  by  the  fact  that  the  propeller  potential  is  more  con¬ 
veniently  expressed  in  cylindrical  coordinates  whereas  the 
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fluid  potential  around  the  ellipsoid  of  revolution  Is  ex¬ 
pressed  In  terms  of  ellipsoidal  coordinates.  This  study, 
therefore,  while  taking  cognizance  of  the  effect  of  the 
boundary  on  the  propeller.  Ignores  the  "image  effect."  It  Is 
believed  that  for  bodies  of  practical  Interest — that  Is,  very 
slender  bodies — the  effect  of  the  so-called  "image  system" 
will  be  secondary.  The  axial  velocity  Induced  by  the  action 
of  the  propeller  Is  much  smaller  than  the  forward  velocity 
of  the  body,  so  that  the  effect  of  the  distortion  In  the 
axial  direction  will  be  small. 

There  are  two  possible  schemes  for  evaluating  the 
axial  vibratory  force  exerted  on  the  prolate  spheroid;  l)  by 
integrating  the  transient  pressure  over  the  surface  of  the 
body}  2)  by  means  of  Laigally*s  theorem  since  the  spheroid  can 
be  represented  by  a  line  source-sink  distribution  of  known 
strength.  It  is  to  be  noted  that  for  a  propeller  located  on 
the  body  axis,  only  the  fore-and-aft  force  is  present;  the 
transverse  force  in  any  direction  will  be  nil  due  to  the 
symmetry. 

According  to  scheme  1,  the  axial  force  AP^  exerted 

on  a  ring  of  radius  r  located  at  x.  will  be  given  by 

o  o 

dx 

AP  =  -  P  sin  a  r  d© 

X  o  cosa 

where  P  is  the  pressure  distribution  arising  from  the  opera¬ 
ting  propeller.  The  resultant  axial  force  exerted  on  the 
ellipsoid  will  be  given  therefore  by 
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by  utilizing  eq.  2.  In  this  form  the  expression  for  the 
pressure  arising  from  the  propeller  Is  written  In  terms  of 
the  nondlraenslonal  coordinate  i  of  the  prolate  spheroid.  The 
pressure  field  emanating  from  the  operating  propeller  is  made 
up  of  a  superposition  of  axially  and  tangentially  directed 
doublets.  The  latter  doublet  distribution  Is  In  a  plane 
normal  to  the  x-axls  and  will  not  contribute  to  the  axial 
component  of  the  vibratory  force,  whereas  the  first  one, 
which  is  designated  as  the  thrust  producing  pressure,  Pip, 
will  be  the  only  contributor  In  eq.  3* 

In  the  second  scheme  the  ellipsoid  of  revolution  Is 
represented  by  a  line  segment  distribution  of  sources  and 
sinks  between  the  focal  points,  of  strength^®  given  by 


where  U  =  forward  velocity  of  the  hull 

Then  the  Interaction  force  between  the  propeller  and  the  ele¬ 
ment  of  the  singularity  representing  the  ellipsoid  will  be 
given  by  application  of  Legally *s  theory  as 

=  -  Wf-*)  (5) 

where 

u  =  axial  component  of  the  velocity  induced  by  the  pro- 
peller  at  any  point  on  the  axis  of  the  ellipsoid. 

It  is  known,  however,  that  the  Instantaneous  linear¬ 
ized  pressure  in  the  field  of  an  operating  propeller  is  given 
by 


P/p 


(6) 


R-855 

-5- 


where  fi  =*  propeller  angular  velocity  and  p  »  fluid  density  and 
0  the  velocity  potential.  The  first  term  is  the  torque- 
associated  pressure  signal;  the  second  is  the  pressure  com¬ 
ponent  associated  with  propeller  thrust.  The  torque-associated 
pressure  acting  in  a  plane  normal  to  the  x-eucls  will  not  con¬ 
tribute  to  the  x-coraponent  of  the  force  as  it  integrates  to 
zero  on  any  circle  with  center  on  the  propeller  axis.  In 
addition,  its  value  on  the  x-axls  is  identically  zero,  since 
r  =0.®  Hence  eq.  6  is  reduced  to 


which  upon  substitution  into  eq.  5  and  integration  leads  to 


f+k 

where  k  =  ae  which  defines  the  locus  of  the  focal  points.  The 
axial  component  of  the  vibratory  force  exerted  on  the  prolate 
spheroid  by  the  propeller  action  is  determined  either  by  eq.  3 
or  eq.  7i  respectively,  once  the  pressure  on  the  hull  or  on 
its  axis  is  determined. 

Pressure  Field 

The  pressure  field  emanating  from  a  marine  propeller 
operating  under  space- varying  inflow  conditions  has  been  the 
subject  of  ref.  3  and  4.  The  flow  about  propeller  blades 
rotating  through  the  space-varying  ship  wake  is  equivalent  to 
that  of  a  wing  moving  through  a  sinusoidal  gust.  Although  the 
problem  is  of  very  pronounced  three-dimensional  character 
(low  aspect  ratio  of  the  blades ,  three-dimensional  gust) ,  the 
two-dimensional  unsteady  aerodynamic  theory  has  been  utilized 
in  a  stripwlse  fashion  as  the  best  expedient  in  the  present 
state  of  the  art.  It  is  assumed  that  a  two-dimensional  flow 
condition  exists  at  every  blade  section  eind  that  three- 


P..  = 


4tt 

TT 


(x-f)  Pn 
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dimensional  effects  can  be  taken  Into  account  by  Burrlll's^^ 
method  which  applies  correction  factors  to  the  axial  and 
tangential  components  of  the  Inflow  velocity  relative  to  the 
propeller. 

At  any  radial  section  of  the  blade  the  velocity 
components  and  the  resultant  velocity  are  shown  in  sketch 
No.  2  below: 


where  V 

n 

a,  a» 
U,  U» 


a* 

a 

<t> 

^o 


propeller  forward  velocity 

8uigular  propeller  velocity 

axial  and  tangential  correction  factors, 
respectively^^ 

resultant  fluid  velocity  with  respect  to  the 
propeller  In  the  absence  and  presence  of  the 
wake,  respectively 

longitudinal  and  tangential  wake  velocity, 
respectively 

angle  of  attack 

effective  angle  of  attack 

hydrodynamic  pitch  angle 

effective  pitch  angle 


R-855 

-7- 


If  the  Instantaneous  velocity  U*  Is  resolved  alo*ig  and  normal 
to  the  chord,  then  the  corresponding  surge  and  cross  veloci¬ 
ties  are  given  as 

-  U  cos  a  - 
s  8 

V-  =  U  sin  a*  +  W 
c  c 

where 

“  a  (V-Wj^)j  cos  T)^  -  jw^  +  a»  (nr-W^)j  sin 

For  small  angles,  the  following  four  combinations 
of  surge  and  cross  velocities  will  simulate  the  flow  around 
the  propeller  and  give  rise  to  corresponding  components  of 
lift  produced  by  a  propeller  operating  In  a  wake: 


Surge  Velocity 

Cross-flow  Velocity 

Lift 

1 

U 

Ua* 

2 

U 

l(2) 

3 

Ua* 

l(3) 

4 

-W 

a 

^c 

The  lift  at  the  corresponding  sections  Is  determined  by 
utilizing  the  results  of  unsteady  aerodynamic  theory  in  terms 
of  the  harmonic  constituents  of  the  wake  flow.  The  elementary 
thrust  and  torque  developed  by  a  propeller  section  at  a  radial 
distance  r^  will  be  given  by 


Trr?  p  I  U«2  sin  X 
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I 


■ 

I 


where 

=  lift  coefficient  per  blade  section 
B  ■=  number  of  blades 

a  =  2^  =  solidity 

b  =  chord  length  at  each  propeller  section 


On  the  assumption  of  constant  blade-loading  distribution 
along  the  chord  the  strength  of  the  pressure  doublet  related 
to  the  thrust  and  torque  will  be  given  by 


^  li 


^  /  .^Bb 
n=  -00 
n/O 


nBb  ^ 
sin  ^  e 


lnBf)(t 


IT 


dT 


(10) 


n=-oo 


nBb 


lnB0(t- 


-L.  ) 

n  anr^' 


1  dQ 

ri  Hri 


The  strengths  of  pressure  doublets  are  finally  expressed  by 
means  of  eqs.9  and  10  in  terms  of  the  lift  coefficients 
associated  with  four  flow  conditions  Indicated  by  eq.  8.  Then 
the  pressure  signals  associated  with  the  thrust  and  torque 
loading  are  given  by  the  well  known  expressions 
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where  r^,6j  =  polar  coordinates  of  the  doublet  singularity  on 
the  propeller  plane 

x,r, ©  =  cylindrical  coordinates  of  the  field  point  or 
point  of  observation 


R  =  propeller  radius 


S  =  Descartes'  distance  between  the  singularity  eind 
field  point  =|x®+r®+ri^  -  2rri  cos(0j^-©)| 


In  ref.  3  the  pressure  signals  emanating  from  the 
propeller  and  associated  with  the  thrust  and  torque  loading 
are  determined  by  means  of  eq.  11  In  conjunction  with 
eqs.  8,  9j  10.  The  thrust  loading  pressure  which  Is  of  in¬ 
terest  In  the  present  work  Is  given  by  (see  eq.  20  for  the 
case  k  =  0  In  the  above  reference) 


.InBfit 


b  r  “  r  “  ^n 

-bJ  C(ri)  S  “2”^  *  n+n  +®3»n-ni)'*‘ 

*'o  L  ni=l  ^ 


Pm  - 


T  471-2 


JL-  b. 


ni=l 


-lx 


(12) 


-^B3,n+ni"®3,n-niJ“  _^(B3,n+ni‘B3,n-ni^ 


dr. 


where 


cos  nB0 


v,n  (x^+r^+r 


pH  _  sin  nB© 
v,n  “  (xa+r^+rf 


-27r 

e_  I  cos  nBg 

I  J  (l-k,COS  5)’'^= 

•'o 

27r 

©  I  cos  nB^  , 

J  (i^.Tossy/= 


(13) 


2r  r 

ki=  -3^ 


,  nBb 


;  c(r,)  =^sln^  e'"^ 


x^+r^+r^  .  --1^  nBb  - 2ri 
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a..  the  cosine  and  sine  Fourier  coefficients  of  the  load- 
n^'  ni 

Ing  function  (eq.  9)  expressed  In  terms  of  the  harmonic  con¬ 
stituents  of  the  wake  velocity.  Equation  13  can  be  written 
In  more  convenient  form  as 


cos  nB0 


_ cos  nBg  dg _ 

[(x®+r®+r®  -2rrj^cos^]^ 


os  nBe  I  cos  nB4  { 


l^^e^cos  m^  fjj^Ck'r) Jjjj(k'ri)a"l*l^  dk’  dS 


0  V  k'=0 

where  use  Is  made  of  the  well  known  Fourier  expansion  of  the 
Inverse  of  the  Descartes*  distance  l/s,  i.e.t 


e„  cos  md  f  J  (k*r)j  (k*r,)e 
tn  I  m  m 


^Ixlk* 


where  =  1  m  =  0 
m 

€„  =  2  for  m  /  0 
m  ' 

Interchanging  the  order  of  summation  auid  Integration  (which 
Is  permissible  In  this  case  as  shown  In  Appendix  2  of  ref.  6) 
leads  to  ,  , 

N  r  -UU' 

xB4,n  =  ”27r  cos  nB0  ^  I  J^(k*r)  J^(k*r  J  e  '  dk’ 

Similarly  „  ,  , 

xB4',n  =  -27rsin  nB0  JnB^  -^nB^ 

Therefore  eq.  12  Is  seen  to  be  made  up  of  terms  of  the  form 

-“olsJn  ^f”  '"‘''''lit’ 

'  •'o 

(eq.  14  con’t  on  next  pg) 
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where  rij^  extends  from  one  to  Infinity. 

However,  the  indicated  ©-integration  in  eq.  3  fixes 
the  order  of  harmonics  which  will  contribute  to  the  blade- 
frequency  vibratory  force.  In  fact,  the  only  terms  which  re¬ 
main  after  the  ©-integrations  are  those  for  n-n^®©  or  ni=n 
since  n/o  and  n+n/O,  both  being  greater  than  zero. 

Therefore,  out  of  ein  infinite  number  of  possible  com¬ 
binations  of  the  harmonics  of  the  space  function  with  those 
of  the  loading  function, only  the  zero  order  harmonic  of  the 
space  function  with  the  blade  frequency  part  of  the  loading 
contributes  to  the  vibratory  force.  Hence  the  blade- frequency 
vibratory  force  referred  to  a  coordinate  system  fixed  at  the 
center  of  the  spheroid  will  be  given  by 


glnEfit  /rR  r  T 

V  '  X  a(l-e2)^^  G(rjj^(a^(rJ-ibj^(ri)J 

00 

J^(k'rj^)j^(k'r^)e'^l^’+^l^'^dk'  dr^d^ 

k=0 

where 

1  /a  i/a 

r;  =  a(l-e2)'^ 

f’  =  f/a 

The  I  -  integration  by  parts  leads  to 


p  =  -  « 

X  “ 


where 


-  7 


Jrmt 


-R 

a(l-e^)  j  C(r,)  [a„(  r,) -lb^(  r,)j 
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f 


where 


or 


or 


[  =  _ 1 _  I  1 _ 

^  C(f'+a)2  +  rf]  Va  [(f-a)®  +  rf] 

■a  =  /*  f  d4 

*^1  •'o 

■  »2/*^  _ I^oJ _ dl  _ _ 

■2  |a2(f+e)=  +  X  (1-^2)  V2  ^  I  £7-2 

1-F)  *^-1/2 


where  K(a)  *=  complete  elliptic  Integral  of  the  first  kind  of 
modulus 


_ 2l/r^X(l-g^)  _ 

{a2(f»+4)2  +  [X(l-e2)  V2  +  ri]2  (  1/2 


X  «  a(l-e®)^/^ 

Q_j.^g  =  Legendre’s  function  of  the  second  kind  of  -I/2 
order 

Z  =  a^(f«+g)^  +  X^(l-g^)  +  r^ 

2rX  (l-?2)^/2 

In  eq.  15  denotes  the  bow  and  stern  contribution  to  the 
vibratory  force  and  gives  the  contribution  of  the  rest  of 
the  surface.  The  second  term  of  gives  the  stern  contribu¬ 
tion,  which  is  the  predominant  contribution,  whereas  the  first 
term  is  associated  with  the  bow  contribution. 

In  the  second  scheme  where  the  pressure  P|j,  has  to 
be  evaluated  at  r  =  0,  eq.  14  indicates  that  0  if  and 

only  if  n  =  0  or  n  ±  n^^  =0.  The  case  n  =  0  is  of  no  Interest 
to  us  since  the  blade- frequency  force  must  be  evaluated.  The 
case  n  +  nj^  =  0  is  impossible  since  n  and  n^  vary  from  1  to 
infinity.  Therefore  the  only  possibility  is  n  -  n^^  »  0  or 
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“  n  which  reduces  eq.  12  to  the  form 


‘  ■  “  ‘CSJI  “  - 

0  •<)  (16) 

After  substitution  of  eq.  l6  into  eq.  7  the  longitudinal  com¬ 
ponent  of  the  vibratory  force  is  determined  as 


^  *^0  4) 


f+ae  R 


p  =  - 

X 


^inBf2t  C  f  , 

V“  - E— j  J  G(rJ  I  a^(rj-lb^(rj|  (f-x) 

*^x=f-ae  0 

/•“ 

^  I  dk'  dr 3^  dx 

which  after  the  k  and  x  integrations  leads  to 


InBfJt  yR 


Px=b: 


'(  f+ae) ^+rf  +  J(f-ae) ®+rf 


f-ae  +w(f-ae)^  +  / 

+  log  ■====z=^  [dr^  (17) 

f+ae  +W(f+ae)^  +  1 

Comparison  of  eq.  15  and  17  shows  that  both  are  of  the  same 
structure  and,  in  fact,  the  first  two  terms  of  eq.  17  can  be 
shown  to  be  identical  with  the  first  two  terms  of  eq.  15  for 
a  slender  body,  with  e— 1  and  ^■^®  ^ .  The  logarithmic 

term  is  equivalent  to  the  term,  which  is  small  and  repre¬ 
sents  the  contribution  of  the  line  singularity  distribution 
extended  between  the  focal  points. 

For  both  expressions,  however,  the  remaining  radial 
integration  must  be  performed  either  numerically,  since 
numerical  values  of  the  loading  function  a^(r3)are  known  in 
terms  of  the  measured  wake  velocity,  or  by  utilizing  the 
mean- value  theorem  of  the  calculus. 
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I 


The  second  scheme,  however,  can  be  further  simpli¬ 
fied  by  assuming  that  the  chordwlse  and  loading  functions  are 
evaluated  at  a  suitable  radial  distance  (r^)  and  then  con¬ 
secutive  Integrations  with  respect  to  x  and  rj^  lead  to  the 
following  form 


_lnBf2t 


Rf  log 


a„(Re)-lb„(R^) 


R(f-fae)  -f  ( f+ae)\/(f-ae)^  +  R^ 
R(f-ae)  +  ( f-ae)  ( f+ae)^  +  R® 


+R^log 


f-ae  +‘ 
f+ae  +■» 


f-ae)^  +  ^ 
f+ae) 2 


The  method  of  evaluating  the  Fourier  coefficients 
a_(r,)and  b  (ri)  has  been  indicated  sketchily  but  the  necessary 
Information  la  presented  In  detail  In  refs,  3  and  4.  With  the 
experience  gained  from  previous  calculation.  It  Is  suggested 


that  the  quasi-steady  approach  for  the  loading  function  be 


utilized.  Making  use  of  the  two-dimensional  unsteady  theory 
In  a  strlpwlse  fashion  has  led  to  a  great  discrepancy  between 
the  measured  and  calculated  vibratory  thrust,  whereas  the 
quasi-steady  results  are  close  to  the  experimental.  This  dis¬ 
crepancy  has  led  to  a  series  of  Investigations  at  Davidson 
Laboratory  which  apply  lifting- surface  theory  approach  to  the 
marine  propeller  case. 

In  the  quasi-steady  flow  case,  then,  the  lift  per 
unit  of  blade  radius  at  each  radial  distance  can  be  written* 


STTpc  K  „  K  (U^a  +  -  Ua  W„  -  W„W  ) 

gs  S  C  8  S  C 


where 


c  «=  semi- chord  at  that  radius 

K  ,K  «  correction  factors  for  cascade  and  finite 
®  ®  aspect  ratio  effects^^ 
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The  remainder  of  the  symbols  are  as  defined  previously. 

It  Is  suggested  that  eqs.  l8  and  19  be  used  In  the 
evaluation  of  the  longitudinal  component  of  the  vibratory 
force.  Numerical  calculations  have  been  obtained  by  means  of 
eq.  l8  to  determine  the  dependence  of  the  longitudinal  vibra¬ 
tory  force  on  propeller  clearance,  as  well  as  on  the  slender¬ 
ness  ratio  of  the  ellipsoid  of  revolution,  for  the  same  load¬ 
ing  conditions  and  a  given  propeller  radius. 


In  the  limiting  case  e-^1,  when  the  ellipsoid  de¬ 
generates  to  a  segment  of  the  x-axls  of  length  2ae  about  the 

origin,  it  is  easily  seen  that  lim  P„  =  0.  When,  however, 

e-1  * 

6-^0  and  the  ellipsoid  degenerates  to  a  sphere  of  radius  a, 
the  vibratory  force  P  Is  obtained  by  three  successive  appll- 

A 

cations  of  L’Hospltal’s  rule  as 


InBfit 


The  space  factor,  {i  +  ^) 


a3R2 

compares  with  -  ,  of  ref.  10  where  the  propeller  action 

(f2+R2)3/a 

is  represented  by  a  sink  disk.  The  additional  term  R^/2f®  Is 
very  small  in  all  cases  of  practical  Interest,  so  that  the 
sink  disk  representation  of  the  propeller  action  is  a  good 
approximation  as  far  as  the  longitudinal  force  is  concerned. 

It  Is  Interesting  to  notice,  however,  that  the  difference 
between  eq.  20  and  that  of  ref.  10  lies  in  the  form  of  the 
loading  distribution.  The  sink  representation  for  the  pro¬ 
peller  action  can  be  proved  exact.  In  the  case  treated 
previously,^®  the  entire  disk  is  activated,  while  In  the 
present  case  only  segments  of  the  disk  participate  in  simulating 
propeller  action.  This  fact  is  the  actual  reason  for  the  ob¬ 
served  differences.  Thus  eq.  20  should  be  considered  more 


R-855 

-16- 


accurate.  In  the  marine  propeller  case  with  very  wide  blades 
the  difference  between  disk  activator  and  activating  sectors 
should  be  small. 

Calculations  Indicate  the  Importance  of  the  pro¬ 
peller  clearance  and  the  slenderness  ratio  of  the  body  and, 
of  course,  the  shape  of  the  afterbody,  since  the  slope  of  the 
cross-sectional  area  curve  depends  on  the  slenderness  of  the 
body. 

In  changing  the  propeller  clearance  from  zero  to 
l/4  of  propeller  diameter  (C  =  .04a),  the  vibratory  force  Is 
reduced  by  about  50^  for  slenderness  ratio  0.1  to  0.2,  which 
Is  a  region  of  practical  Interest.  Figure  1  Is  a  chart  of 
the  space  function  (of  eqs.  l8  and  20)  versus  slenderness 
ratio  ^  =^1  -  e^  for  R  =  0.08a. 

CONCLUSION 

In  the  foregoing  analysis  an  expression  for  the 
longitudinal  component  of  the  vibratory  force  exerted  on  an 
ellipsoid  of  revolution  by  the  propeller  action  has  been 
developed  In  closed  form  in  terms  of  elementary  functions,  and 
the  relative  Importance  of  the  various  parameters  has  been 
revealed.  It  Is  Indicated  that  the  vibratory  force  sorongly 
depends  on  the  axial  propeller  clearance  as  well  as  on  the 
slenderness  of  the  ellipsoid  of  revolution  and  therefore  on 
the  slope  of  the  cross-sectional  area  at  the  afterbody  section. 
In  fact.  In  the  region  of  practical  Interest,  for  slenderness 
ratio  of  0.1  to  0.2,  changing  the  propeller  cleareince  from 
zero  to  l/4  propeller  diameter  reduces  the  longitudinal  vibra¬ 
tory  force  to  approximately  half.  In  this  range,  also,  the 
force  decreases  rapidly  with  increasing  fineness. 
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FIGURE  i.  VARIATION  OF  SFVVCE  FUNCTION  WITH  SLENDERNESS  OF  SPHEROID  AT  VARIOUS 
PROPELLER  CLEARANCES 
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